Elliptic Flow Analysis at RHIC with the Lee- Yang Zeroes Method in a Relativistic 

Transport Approach 



O 
O 

(N 

C 

a 



> 

ON 

o 
o 
o 



o 

=5 



X 



Xianglei Zhu, 1,2,3 Marcus Bleicher, 2 and Horst Stocker 1,2 

'Frankfurt Institute for Advanced Studies (FIAS), 
Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany 
2 Institut fur Theoretische Physik, Johann Wolfgang Goethe- Universitat, 
Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany 
'Physics Department, Tsinghua University, Beijing 100084, China 

The Lee- Yang zeroes method is applied to study elliptic flow (1)2) in Au+Au collisions at ^/s — 
200A GeV, with the UrQMD model. In this transport approach, the true event plane is known 
and both the nonflow effects and event-by-event V2 fluctuations exist. Although the low resolutions 
prohibit the application of the method for most central and peripheral collisions, the integral and 
differential elliptic flow from the Lee- Yang zeroes method agrees with the exact V2 values very well 
for semi-central collisions. 

PACS numbers: 25.75.Ld, 25.75.Dw, 25.75.Gz 



Anisotropic flow 0, |E IE 13 > ° r more specifically, ellip- 
tic flow (^2), which is the second Fourier harmonic |5| in 
the transverse distribution of the emitted particles, is ex- 
pected to be sensitive to the early pressure gradients and 
therefore to the equation of state Q50S) of the formed 
fireball in heavy- ion collisions [E IE IE IE IE E3 ■ In addi- 
tion, the elliptic flow of high px particles is related to jet 
fragmentation and energy loss of the primordially pro- 
duced hard antiquark-quark pair when traveling through 
the hot QCD medium |ll|. Therefore, understanding the 
microscopic nature of the elliptic flow is of great impor- 
tance for the study of the key questions of experimental 
and theoretical relativistic heavy ion collisions, the trans- 
port coefficients and the pressure, namely the equation 
of state. 

To measure flow, experiments usually use the reaction 
plane method fl^lllp o r the equivalent two-particle cor- 
relation method~ |l2l Il4 Ho| . However, these two-particle 
correlation-based methods suffer from correlations un- 
related to the reaction plane. These additional contri- 
butions are usually dubbed nonflow effects [ly, such as 
the overall transverse momentum conservation, small an- 
gle azimuthal correlations due to final state interactions, 
resonance decays, jet production 1 171 and quantum cor- 
relations due to the HBT effect [lg- I n order to de- 
crease the contribution of the nonflow effects to the flow 
measurements, many-particle cumulant method was pro- 
posed 0]. In this method, usually 4- and 6-particle cu- 
mulants are used to estimate the collective flow. This 
method has been a ppl ied to the flow analysis at RHIC. 

sps and sis/gsi EE E3, H3, H3, 13 H3, 

And it has also been tested in the UrQMD model in [27j • 
It was found that the two-particle correlation method is 
strongly affected by nonflow effects in the UrQMD model, 
such as jet production and resonance decays, while the 
4- and 6-particle cumulants give the true v 2 values with 
better than < 5% accuracy for the semi-central collisions 
at RHIC, even though this method is affected severely by 
the rather large vi fluctuations for the central and very 
peripheral collisions . More recently, the Lee- Yang ze- 



roes method was developed to analyze the collective flow 
Ull Llii Liil • This new method enables the extraction of 
the "true" collective flow from the genuine correlation 
between a large number of particles. It is expected to 
give the cleanest values of the genuine collective flow. It 
is, in fact, also easier to apply to data than the cumulant 
method. This method has been applied in the flow anal- 
ysis at SIS/GSI 26j. Here it was found that the V\ and 
i>2 values extracted from this new method do agree with 
those from the higher order cumulant method. 

In the present paper, the Lee- Yang zeroes method will 
be tested in the integral and differential v 2 analysis of 
UrQMD events at RHIC energy. The UrQMD model 
32,123 is well suited for the test: it has been shown in 
27 1 that both the nonflow effects and the v 2 fluctuations 
, which can affect the flow data, do exist naturally 
in the model. Most importantly, the true reaction plane 
angle is known by construction in the UrQMD model. 
This allows for the direct calculation of the exact elliptic 
flow from its basic definition, i.e. v 2 — (cos2(cf> — 
Through the present analysis, it will be demonstrated 
that the Lee- Yang zeroes method works well for the v 2 
measurements at RHIC energy. 

In order to determine the "integrated" elliptic flow, 
which is defined as 

V 2 = (^J^i - $r))J , 

the generating function of the Lee- Yang zeroes method 
is introduced and defined as 129, 13CI l31l 



G e {ir) = / JJ [1 + iruj cos ( 



0)} 



In the above definitions, M is the number of particles 
used to estimate the integral flow in each event, <f>j is 
the particle azimuthal angle, r is a positive real variable, 
9 is an arbitrary reference angle and uij is the particle 
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weight. It is worth stressing that in our analysis, the 1/M 
weight is used to decrease the effect of event- by- event M 
fluctuations. Therefore, V2 — «2 in our analysis. Fig. 
shows the amplitudes of the generating function |G e (ir)| 
as a function of r for 9 — at different centralities. For 
the semi-central bins (e.g. 20-30%), with the increase 
of r, \G e (ir)\ decreases rapidly from the value of 1 at 
r = 0, then reaches a sharp minimum which is very close 
to 0. This appearence of the sharp minimum indicates 
that there exists a zero of the generating function. After 
finding the position r$ of the first zero of \G e (ir)\, the 
integral elliptic flow is determined as 

«f =joi/rg, 

where joi = 2.40483 is the first root of the Bessel function 
Jo(x). The final value of the integral V2 is the average 
of v\ over 5 equally spaced values of 9 from to 47r/5, 
and the statistical uncertainty of the final integral vi is 
about a factor of 2 smaller |3C| . 
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FIG. 2: (Color online) The UrQMD model results for integral 
V2 values for Au+Au collisions at y/s — 200A GeV. Results 
from the Lee- Yang zeroes method (v2{oo}) with unit weight 
and 1/M weight are compared to the exact «2 in different 
centrality bins, (see text for details.) 
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FIG. 1: (Color online) \G e (ir)\ as a function of r for 9 = at 
different centralities. 

Fig. |2 shows the Lee- Yang zeroes results on the cen- 
trality dependence of the integral v%. For the integral 
V2 analysis, all particles in the pseudorapidity region 
\r)\ < 2.5 are used and the number of particles from event 
to event fluctuates in each centrality bin. The centralities 
in our analysis are selected according to the same geomet- 
rical fractions of the total cross section (0-5%, 5-10%, 10- 
20%,20-30%,30-40%,40-50%,50-60%,60-70%) as used by 
the STAR experiment jy). Here, however, we use im- 
pact parameter cuts instead of multiplicity cuts. More 
than 1.3 • 10 6 minimum bias events are used in the inte- 
gral V2 analysis. 

The V2 values extracted from the Lee- Yang zeroes 
method with 1/M weight (solid circles) agrees with the 
exact V2 (open triangles) very well for the semi-central 
bins (about 5-40%), as has been seen in Fig. [21 How- 
ever, the Lee- Yang zeroes method does not work well 
for the most central (0-5%) and peripheral (>40%) bins 
in the present model study: the resolution parameter x, 
which is related to the event-plane resolution (3j| , decides 
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FIG. 3: The UrQMD model results for «2(pr) in semi-central 
(20-30%) Au+Au collisions at y/1 = 200A GeV. Results from 
the Lee- Yang zeroes method (t>2{oo}) are compared to the 
exact V2- 

whether the method can be applied or not [H HJ. If 
this parameter is too small, i.e., x < 0.5, the statistical 
error will be too large, and the method is inapplicable, x 
is roughly given by i^v'M. Therefore, it will be smaller 
for the most central bin where V2 is small, and for pe- 
ripheral bins where M is small. For the feasibility test 
with the UrQMD model, x is about 1 in the semi-central 
bin (20-30%), while it is close to 0.5 in the most central 
(0-5%) and peripheral bins (>40%). In fact, as is shown 
in Fig. ^ the first minimum of \G {ir)\ is not compat- 
ible with a zero for most central (0-5%) and peripheral 
(>40%) bins. Therefore, the Lee- Yang zeroes method 
can not be used in most central (0-5%) and peripheral 
(>40%) bins from the UrQMD model. 

We have also tried the Lee- Yang zeroes method with 
unit weight. These results are shown as open circles in 



3 



Fig. El For the semi-central bins (10-40%), it slightly 
undershoots the exact v 2 , which indicates that the ze- 
roes method with unit weight is not completely free from 
the effect of M fluctuations. However, the effect of M 
fluctuations is small. Experimentally, the centrality is 
usually selected according to the multiplicity M of the 
events. In this case, the M fluctuations will be narrower, 
and the effect of M fluctuations on the zeroes method 
will be weaker. The present Lee- Yang zeroes analysis 
based on the M selected centralities (not shown in this 
Letter) indeed tells that the integral i> 2 values from the 
1/M weight and unit weight agree well with each other 
and the M fluctuations have almost negligible effects. 
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trality bin with better than 2% accuracy, However it is 
still necessary to see whether it reproduces the differen- 
tial V2 correctly. Especially at large transverse momenta 
(pr), nonflow contributions are expected to be large. Fig. 
[3] and Fig. 0] shows the comparisons of the Lee- Yang ze- 
roes results for V2{pt) and V2{ff) to the exact i> 2 values. It 
is clear that the results from the Lee- Yang zeroes method 
agree with the exact i> 2 values very well for the whole 
pr- and r;-range. This shows that the Lee- Yang zeroes 
method indeed reproduced exactly the genuine elliptic 
flow values for these centrality bins very well. 

To summarize, we have applied the Lee- Yang zeroes 
method in the V2 analysis of the UrQMD model at 
200 A GeV RHIC Au+Au collisions. Only for the semi- 
central bins (5-40%), is this method actually applicable. 
Both the integral and differential v 2 from the Lee- Yang 
zeroes method agree well with the exact values in these 
bins. The UrQMD model includes many of the nonflow 
correlations which are also seen in the real data. Hence, 
these good agreements show that the Lee- Yang zeroes 
method is not sensitive to the nonflow correlations and 
V2 fluctuations for the semi-central bins and indeed gives 
the genuine elliptic flow there. It should be noted that the 
current UrQMD model underpredict the RHIC v 2 data 
by about 40%. Therefore, in the experiments at RHIC 
or LHC, the zeroes method will work over a somewhat 
broader range of centrality bins. 



FIG. 4: The UrQMD model results for V2(rj) in semi-central 
(20-30%) Au+Au collisions at </a = 200A GeV. Results from 
the Lee- Yang zeroes method (i>2{oo}) are compared to the 
exact V2. 

After determining 7p, the differential w 2 can be esti- 
mated from Eq. (9) in [21| . For the test of the differential 
flow, we use more than 6 • 10 5 semi-central events (with 
impact parameters from 6.7 to 8.3 fm, corresponding to 
about 20% to 30% of the total cross section). From the 
above results on the integral v 2 , we know that the Lee- 
Yang zeroes method reproduces the exact v 2 in this cen- 
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